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ABSTRACT 



Context. Ragozzine & Brown [2007] presented a list of candidate members of the first collisional family to be found among the 
trans-Neptunian Objects (TNOs), the one associated with (136108) Haumea (2003 EL6i). 

Aims. We aim to identify which of the candidate members of the Haumea collisional family are true members, by searching for water 
ice on their surfaces. We also attempt to test the theory that the family members are made of almost pure water ice by using optical 
light-curves to constrain their densities. 

Methods. We use optical and near-infrared photometry to identify water ice, in particular using the (J - Hs) colour as a sensitive 
measure of the absorption feature at 1.6 \im. We use the CH 4 filter of the new Hawk-I instrument at the VLT as a short //-band (H s ) 
for this as it is more sensitive to the water ice feature than the usual H filter. 

Results. We report colours for 22 candidate family members, including NIR colours for 15. We confirm that 2003 SQ317 and 2005 
CB79 are family members, bringing the total number of confirmed family members to 10. We reject 8 candidates as having no water ice 
absorption based on our Hawk-I measurements, and 5 more based on their optical colours. The combination of the large proportion of 
rejected candidates and time lost to weather prevent us from putting strong constraints on the density of the family members based on 
the light-curves obtained so far; we can still say that none of the family members (except Haumea) require a large density to explain 
their light-curve. 

Key words. Kuiper Belt; Methods: observational; Techniques: photometric; Infrared: solar system 



1. Introduction 

The trans-Neptunian object (TNO) (136108) Haumea (2003 
EL 61 ) was discovered by Santos-Sanz et al. [2005] and quickly 
attracted a lot of attention as a highly unusual body. It is one 
of the largest TNOs [Rabinowitz et al., 2006; Stansberry et al., 
2008] and yet is a fast rotator (period ~ 3.9 hours) with a highly 
elongated shape [Rabinowitz et al., 2006]. Its surface was shown 
to be dominated by water ice by Near Infra-Red (NIR) spec- 
troscopy [Tegler et al., 2007; Trujillo et al., 2007; Merlin et al., 
2007; Pinilla-Alonso et al., 2009], yet has a high density of 
2.5-3.3 g cm -3 [Rabinowitz et al., 2006]. It was found to have 
two satellites [Brown etal., 2005, 2006], which also have wa- 
ter ice surfaces [Barkume et al., 2006; Fraser & Brown, 2009]. 
Lacerda et al. [2008] found that Haumea presents hemispherical 
colour heterogeneity, with a dark red 'spot' on one side, using 
high precision photometry. 

Brown et al. [2006] and Barkume et al. [2006] postulated 
that the density, shape and water ice surface could be explained 
by a large collision early in the history of the Solar System. 
Brown etal. [2007] then identified a family of 6 TNOs (1995 
SM55, 1996 TO 66 , 2002 TX 300 , 2003 OP 32 and 2005 RR43), in 
addition to Haumea and its satellites, with orbits that could be 
linked to Haumea and water ice surfaces, which were also at- 
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tributed to coming from this massive collision. This theory re- 
quired that the proto-Haumea was a very large body (radius 
~ 830 km) that had already differentiated early in the formation 
of the Solar System, and that the collision stripped nearly all of 
the outer (water ice) mantle (~ 20% of the total mass of the orig- 
inal body). This left the dense core as Haumea with a thin coat- 
ing of water ice and created a family of re-accumulated lumps of 
almost pure water ice. Ragozzine & Brown [2007] find that the 
collision must have taken place in the early Solar System (with 
an age of at least 1 Gyr), although the lack of weathering on the 
surfaces may imply young bodies [Rabinowitz et al., 2008]. The 
existence of such a family has implications for the dynamics of 
the Kuiper Belt [Levison et al., 2008]. 

Ragozzine & Brown [2007] performed a dynamical study 
and identified two further family members (2003 UZin and 
1999 OY3) with strong dynamical links to the family and colours 
consistent with water ice, and also published a list of candi- 
date family members that had orbital elements consistent with 
this dynamical family, totalling 35 objects including the known 
members. Most of these candidates lacked the NIR spectra 
that could identify water ice on their surfaces though, so they 
remained only potential family members. The diffusion time 
and interaction with resonances make it possible for interlop- 
ers to appear close to the family dynamically, so it is essen- 
tial to have both dynamical and physical properties characterisa- 
tion to confirm family membership [Cellino et al., 2002]. Some 
could be ruled out by either existing NIR spectra (Makemake 
has a methane ice surface; Dumas et al. [2007]; Brown et al. 
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[2007]) or by very red optical colours (1996 RQ 20 , 1999 CD H8 , 
1999 KR 16 , 2002 AW 197 , 2002 GH 32 ; see table 4 for refer- 
ences) or a strong red slope in optical spectra (2005 UQ513; 
Pinilla-Alonso et al. [2008]). Schaller & Brown [2008] subse- 
quently published NIR spectra which confirmed 2003 UZni and 
2005 CB79 as family members, and rejected 2004 SB60- We ob- 
served 13 of the 18 remaining candidate objects (along with 
some of the already characterised objects) with the goal of pro- 
viding this physical information, to identify those with water ice 
surfaces and also to test the idea that these family members could 
be made of nearly pure water ice. We describe our observations, 
the results from them, and their implications in the following 
sections. 

2. Observations and data reduction 

The best method to test for water ice on the surface of a Solar 
System body is through NIR spectroscopy, as water ice has 
strong absorption bands at ~ 1.6 and ~ 2.0 ftm, but this is 
only possible for the brightest TNOs (K < 18). Still, it is possi- 
ble to get an indication of the presence or absence of water ice 
for fainter bodies using photometry, which can be performed on 
smaller (fainter) TNOs. 

We conducted the observations at the European Southern 
Observatory (program IDs: 81.C-0544 & 82.C-0306), on both 
the La Silla and Paranal (VLT) sites. Observations in the visible 
wavelengths (BVRi filters) were performed using the EFOSC2 
instrument [Buzzoni et al., 1984] mounted on the NTT [since 
April 2008; Snodgrass et al., 2008]. This is a focal reducing im- 
ager and spectrograph with a single CCD. The near-infrared 
observations (J, CH4 bands) were performed using the newly 
commissioned wide-field camera Hawk-I [Pirard et al., 2004; 
Casali et al., 2006]. We had three observing runs scheduled with 
each instrument, as detailed in table 1. This table lists all ob- 
jects we attempted to observe, although not all were detected and 
some time was lost to poor weather conditions. In particular the 
June 17th Hawk-I run (run B) was very badly affected by clouds, 
with only 1999 KRie reliably detected in both bands. Exposure 
times were generally 300-600 seconds in the optical, while in the 
NIR we took sequences of J-CH4-J to give an average J magni- 
tude at the time of the CH4 observations, and to confirm identifi- 
cation of the object based on its motion between the two sets of 
/-band images. The CH4 filter observations took the largest part 
of the time; between 15 minutes for the brightest objects to a few 
hours for the faintest ones, each split into short individual expo- 
sures and dithered due to the bright NIR sky. Note that due to 
the long effective exposure times any variation (due to shape or 
albedo variation across the surface) is smeared out, and cannot 
be detected in our NIR data. 

The advantage of using Hawk-I is that the CH4 band filter 
is a medium width filter with a wavelength range that is entirely 
within the broad water ice absorption between 1 .4 and 1 .75 fim. 
The standard //-band is broader and covers a range that is part in 
and part out of this band 1 . We therefore use the CH4 filter as a 
short H filter (henceforth H$ ) which gives a colour measurement 
(J - Hs) that is very sensitive to water ice absorption. All of the 
filters used in this work are listed in table 2. 

The data were reduced in the normal manner (bias sub- 
traction, flat fielding, sky subtraction etc. as appropriate). For 
the EFOSC2 data the objects were generally visible in individ- 
ual frames and aperture photometry was performed directly on 

1 See http://www.eso.org/sci/facilities/paranal/instruments/hawki/inst 
for transmission curves. 



Table 1. Observational circumstances. 



Object 7 /Y 5 a 7 Run 3 Epochs 7 
(#) (Designation) (AU) (AU) (°) BVRi 



48.8 48.0 0.7 D 



Notes: 

" Heliocentric distance; b Geocentric distance; c Phase angle. 
d Runs: A = 2008 June 3rd - 5th, EFOSC2; B = 2008 June 17th, 
Hawk-I; C = 2008 August 30th - September 1st, EFOSC2; D = 2008 
September 9th, Hawk-I; E = 2008 December 29th - 31st, EFOSC2; F 
= 2009 January 4th, Hawk-I. 

' Number of epochs observed in each filter (for EFOSC2 runs). 



each, using the optimum aperture based on the measured stellar 
FWHM in each frame and an average aperture correction mea- 
sured using the field stars [see Snodgrass et al., 2005]. Where 
multiple epochs were obtained we then report a weighted mean 



20161 
181855 
40314 

86047 
86177 
130391 

55565 



55636 
136108 



120178 



136199 



120347 



202421 



1996 RQ 20 


39.6 


39.0 


1.1 


C 






4 




1996 TR 66 


40.3 


40.0 


1.4 


E 


2 


2 


2 


2 


1998 HL 151 


38.9 


38.2 


1.0 


A 


2 


2 


2 


2 


1998 WT 31 


38.0 


37.3 


1.0 


E 


2 


2 


10 


2 


1999 CD 158 


47.6 


46.5 


0.6 


E 


1 


1 


24 


1 


1999 KR 16 


36.3 


35.6 


1.2 


B 










1999 OH 4 


39.1 


39.6 


1.3 


A 


1 


1 


1 


1 


» 


39.1 


38.2 


0.6 


C 


2 


2 


2 


2 


1999 OK 4 


46.4 


45.8 


1.1 


A 


1 


1 


1 


1 


1999 OY 3 


40.1 


39.7 


1.3 


A 






11 




» 


40.2 


39.5 


1.1 


B 










» 


40.2 


39.4 


0.8 


D 










1999 RY 215 


35.8 


34.8 


0.2 


C 






21 




» 


35.8 


34.8 


0.3 


D 










2000 CG 105 


46.8 


46.1 


0.8 


E 


2 


2 


22 


2 


2000 JG 81 


34.8 


33.8 


0.5 


A 


1 


1 


1 


1 


2001 FU 172 


31.8 


30.9 


1.0 


A 


1 


1 


1 


1 


2001 QC 298 


40.6 


39.6 


0.3 


C 






17 




» 


40.6 


39.6 


0.2 


D 










2002 AW 197 


46.6 


45.8 


0.7 


F 










2002 GH 32 


43.1 


42.2 


0.7 


A 


3 


3 


18 


3 


» 


43.1 


42.4 


1.0 


B 










2002 TX 300 


41.4 


40.6 


0.8 


D 










Haumea 


51.1 


50.6 


1.0 


A 


1 


1 


1 


1 


» 


51.1 


50.8 


1.1 


B 










» 


51.1 


51.1 


1.1 


F 










2003 HA 57 


32.7 


32.0 


1.3 


A 


1 


1 


1 


1 


» 


32.7 


32.2 


1.6 


B 










2003 HX 56 


46.5 


45.9 


1.0 


A 


2 


2 


2 


2 


2003 OP 32 


41.4 


40.6 


0.6 


D 










2003 QX91 


33.6 


32.6 


0.5 


C 






4 




2003 SQ317 


39.3 


38.3 


0.6 


C 






15 




» 


39.3 


38.3 


0.4 


D 










2003 TH 58 


36.0 


35.1 


0.5 


E 


2 


2 


23 


2 


» 


36.0 


35.1 


0.7 


F 










Eris 


96.7 


95.9 


0.4 


D 










2003 UZ 117 


39.4 


38.9 


1.3 


D 










2004 PT 107 


38.3 


37.9 


1.4 


A 


4 


4 


24 


4 




38.3 


37.7 


1.3 


B 












38.3 


37.4 


0.7 


D 










2004 SB 60 


44.0 


43.1 


0.6 


C 






16 




2005 CB 79 


40.1 


39.3 


0.9 


E 


1 


1 


2 


1 




40.0 


39.2 


0.8 


F 










2005 GE 187 


30.8 


29.9 


0.9 


A 


3 


3 


33 


3 




30.8 


30.1 


1.3 


B 












30.8 


31.1 


1.7 


C 






17 






30.8 


31.3 


1.6 


D 










2005 UQ513 


48.8 


48.1 


0.8 


C 
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Table 2. Filters used in this study. 



Filter 


Instrument 


A c 


AA 






yum 


yum 


B 


EFOSC2 


0.440 


0.094 


V 


EFOSC2 


0.548 


0.113 


R 


EFOSC2 


0.643 


0.165 


i 


EFOSC2 


0.793 


0.126 


J 


Hawk-I 


1.258 


0.154 


B. s (CH4) 


Hawk-I 


1.575 


0.112 



Notes: 

A c = Central Wavelength, AA - 



Bandwidth. 
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We used this relation to generate the expected colour, and there- 
fore H s magnitude, for each 2MASS star in each field, which 
were then used to give the calibrated Hs magnitude for the 
TNOs. We also used the same approach to derive the colour 
term for the difference between 2MASS and Hawk-I J bands, 
and found that the Hawk-I J does not significantly differ from 
the 2MASS band, as expected. We note that the spectral types 
further from the linear trend fall into two groups; those below 
the trend at (J - #)2M * are B stars that do not feature in our 
NIR images, while the 'tail' that curves away from the line at the 
red end is made up of M giants, with M8-10 being significant 
away from the linear relation. These are separable from the rest 
of the sample though as giants have a very red 2MASS (J - K) 
colour; Brown [2003] show that stars with (J-K) > 0.5 are most 
likely giants, while we find that using limit of (J — K) < 1 .26 re- 
moves the M8-10iii stars that do not fit the linear trend while 
keeping other stars. Having said this, we note that the exclu- 
sion or inclusion of these stars made no significant difference to 
our calibration as there were very few late M giant stars within 
our sample. The colours of the 2MASS stars in the fields ob- 
served were approximately normally distributed around a mean 
(J - H)j\k - 0.6 with a standard deviation of 0.2. 



ro - 



3. Colours 



^ ^2MASS 

Fig. 1. Theoretical difference between 2MASS H and Hawk-I 
H s for different stellar spectra, as a function of 2MASS (J - H). 



magnitude. This approach allowed us to look for variation in the 
/?-band magnitude for those objects where we obtained a light- 
curve. For fainter objects the images were shifted based on the 
predicted motion of each object and combined to give a deep im- 
age per filter. We also produced equivalent combined images of 
the star fields (no shifts) in which we could measure the bright- 
ness of field stars for photometric calibration. For Hawk-I all 
data were shifted and combined as the individual exposures were 
short because of the high sky background in the NIR. 

The EFOSC2 data were calibrated in the normal way, via 
observations of standard stars from the Landolt [1992] cata- 
logue. The EFOSC2 /-band data was calibrated directly onto the 
Landolt scale; this filter is very close to the standard Cousins I- 
band used by Landolt. Data from non-photometric nights were 
calibrated via observation of the same fields on later photometric 
nights, to calibrate the field stars as secondary standard stars. 

Calibration of the Hawk-I data was a more involved process 
as it contained the non-standard filter Hs- The J and H band 
magnitudes of all available stars in each field were taken from 
the 2MASS point source catalogue [Skrutskie et al., 2006]. We 
then generated theoretical colours (H 2 m - Hs) for stars of all 
spectral types (O-M) by convolving the response of the 2MASS 
H and the Hawk-I Hs with spectra from the libraries of Pickles 
[1998] and Ivanovetal. [2004] : . For stars the resulting differ- 
ence is linearly related to the 2MASS ( J - H) colour (Fig. I ): 

(H 2M - H s ) = -0.097(7 - H) 2M - 0.019. (1) 

2 These libraries can be downloaded from the ESO web pages at 
http://www.eso.org/sci/observing/tools/standards/IR_spectral_library_new/ 



We report the resulting photometry in table 3, where we give 
the mean magnitude in each band at each epoch and also an in- 
dication of the variation seen in the /?-band where we obtained 
light-curves. In table 4 we give the average colours of all fam- 
ily members that have published photometry, including our own 
results, taking a weighted mean where multiple measurements 
exist. From these average colours we calculate reflectances by 
comparing them to the Solar colours. To calculate the reflectance 
in the Hs band we used a theoretical (J - Hs) colour for the 
Sun generated by convolving the response of these filters with 
the Solar spectrum. We subsequently confirmed this value by 
observing a Solar analogue star with Hawk-I: The theoretical 
(J - Hs) & = 0.273, while the value measured for the Solar 
twin S966 (taken from the catalogue of Solar twins in M67 by 
Pasquini et al. [2008]) is (J - H S ) Q = 0.288 + 0.007. These 
are consistent at the level of the uncertainty on our TNO colour 
measurements. We also report the visible slope for each object 
(%/100 nm) in table 4, calculated from the reflectances via a lin- 
ear regression over the full BVRI range when it is available, or 
whichever measurements exist in other cases. 

The reflectance 'spectra' of the TNOs from this photometry 
are shown in Fig. 2, for all objects with photometry in at least 
three bands. The combined visible and NIR spectrum of Haumea 
from Pinilla-Alonso et al. [2009] is shown for comparison to the 
photometry. The large TNOs Eris (not a family member; ob- 
served for comparison) and Makemake (dynamically a family 
member candidate) are known to have methane ice surfaces from 
NIR spectroscopy [Dumas et al., 2007; Brown et al., 2007] and 
clearly differ from the Haumea spectrum. Note that those ob- 
jects marked with an asterisk in the figure have their reflectance 
normalised to the R -band, as no V-band photometry was avail- 
able. For Haumea-like neutral spectra this makes no difference, 
but this could give an offset in the case of red slopes; these four 
spectra should not be directly compared with the others in the 
figure, but can be compared with the Haumea spectrum. 
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Table 3. Photometry. Mean apparent magnitudes for each object at each epoch. 





Object 


Run" 




B 


V 


R 








J 


Hs 


Am/ 




1996 RQ 20 


C 




- 


- 


22.95 ± 0.05 




- 




- 


- 


- 




1998 HL 151 


A 


25 


37 ± 0.28 


24.25 ±0.12 


23.87 ±0.13 


23 


.08 ± 


0.17 


- 


- 


- 


181855 


1998 WT 31 


E 


24 


.52 + 0.15 


23.81 ± 0.11 


23.24 ± 0.06 


22 


.69 ± 


0.14 


- 


- 


<0.1 




1999 CD 158 


E 


23 


.08 ± 0.07 


22.31 ±0.06 


21.68 ± 0.01 


21 


.20 ± 


0.06 


- 




0.6 


40314 


1999 KR 16 


B 






- 


- 




- 




20.02 ± 0.07 


19.47 ± 0.10 


- 




1999 OH 4 


C 


25 


.01 ± 0.22 


22.39 ± 0.09 


22.19 ± 0.10 


21 


.76 ± 


0.32 


- 


- 


- 


86047 


1999 OY 3 


A,D 




- 


- 


22.26 ± 0.03 




- 




21.78 ± 0.10 


22.04 ± 0.35 


<0.1 


86177 


1999 RY 215 


C,D 




- 


- 


22.16 ±0.01 




- 




21.19 ± 0.14 


20.66 ± 0.17 


<0.1 




2000 CG 105 


E 


24 


.14 + 0.09 


22.60 ± 0.04 


22.62 ± 0.02 


22 


.52 ± 


0.09 






0.45 




2001 FU 172 


A 


26 


.71 ± 1.52 


23.80 ±0.15 


23.13 ± 0.12 


22 


.66 ± 


0.19 


- 


- 


- 




2001 QC 298 


C,D 






- 


22.18 ± 0.03 




- 




21.16 ±0.08 


20.65 ±0.12 


0.4 


55565 


2002 AW 197 


F 




- 


- 


- 




- 




18.50 ± 0.05 


18.11 ±0.06 


- 




2002 GH 32 


A 


23 


.91 ±0.09 


21.87 ±0.05 


21.87 ± 0.02 


19 


.96 ± 


0.09 


- 


- 


0.75 


55636 


2002 TX 300 


D 






- 


- 




- 




18.67 ± 0.07 


19.14 ± 0.10 


- 


136108 


Haumea 


F 




- 


- 


- 




- 




16.46 ± 0.07 


17.06 ± 0.08 


- 




2003 HX 56 


A 


25 


.25 ± 0.36 


24.03 ±0.16 


23.68 ± 0.16 


23 


.42 ± 


0.43 


- 


- 


- 


120178 


2003 OP 32 


D 




- 


- 


- 




- 




19.08 ± 0.05 


19.58 ± 0.06 


- 




2003 QX91 


C 




- 


- 


23.66 ±0.12 




- 






- 


- 




2003 SQ317 


C,D 




- 


- 


22.05 ± 0.02 




- 




21.59 ± 0.05 


22.04 ± 0.19 


1.0 




2003 TH 58 


E 


23 


.50 ± 0.05 


22.89 ± 0.04 


22.51 ± 0.02 


22 


.03 ± 


0.04 


21.73 ±0.09 


20.45 ±0.18 


<0.1 




2003 UZ 117 


D 
















20.24 ± 0.07 


20.86 ±0.10 






2004 PT 107 


A,D 








21.66 ± 0.01 








20.41 ± 0.14 


19.87 ±0.18 


0.05 


120347 


2004 SB 60 


C 








20.21 ± 0.01 












0.2 




2005 CB 79 


E,F 


21 


.45 ± 0.02 


20.71 ± 0.03 


20.36 ± 0.02 


19 


.98 ± 


0.03 


19.67 ± 0.07 


20.18 ± 0.16 






2005 GE 187 


A 


23 


.76 + 0.10 


22.78 ± 0.09 


22.02 ± 0.01 


21 


.47 ± 


0.11 






<0.1 






C,D 








22.13 ±0.03 








20.84 ± 0.08 


20.18 ± 0.12 


0.5 


202421 


2005 UQ513 


C,D 








20.30 ± 0.01 








18.89 ± 0.07 


18.59 ± 0.10 


0.3 


136199 


Eris 


D 
















17.73 ± 0.07 


17.49 ± 0.09 





Notes: 

" Runs A-F as listed in table 1 . 

* Am R is the variation in i?-band magnitude seen for objects where (partial) light-curves were obtained. The uncertainty on each is ~ 0.1 mag. 



4. Discussion 

4.1. Family membership 

We first wish to determine which candidates are actually fam- 
ily members, and which are dynamical interlopers with different 
surface properties. We find that the (J - Hs) colour is a good 
diagnostic of the presence or absence of the water ice absorp- 
tion feature at 1.6 fim, as expected: For Haumea we measure 
(J - Hs) = -0.60 ± 0.11, and the colour is also significantly 
negative for the other known family members observed, while 
for the methane ice dominated comparison TNO Eris we find 
(J - Hs) = 0.25 ±0.11. The colours for all objects are given in 
table 4, along with the visible slopes, and these are also plotted in 
Fig. 3. In the figure there is a clear separation between the fam- 
ily members with negative (J - Hs) at the bottom and the other 
objects at the top, and also a tendency for those with water ice to 
have blue/neutral surfaces (shallower slopes). While those with- 
out water ice have a large range of slopes from neutral to very 
red, there are no bodies in the lower right of the figure (water ice 
and red slope). We use this separation to make a rough assess- 
ment of the family membership for candidates with only optical 
colours; we can rule out membership for objects with very red 
slopes, but cannot use a blue slope to confirm membership. 

We confirm two more family members in addition to those 
listed by Ragozzine & Brown [2007]; 2003 SQ317 and 2005 
CB79. These have (J - H s ) = -0.45 + 0.20 and -0.50 ±0.17 
respectively. 2005 CB79 has since been confirmed as a family 
member by NIR spectroscopy [Schaller & Brown, 2008]. For 
2003 SQ317 the lack of optical colours as supporting evidence 



and the relatively large uncertainty on (J - Hs) makes the wa- 
ter ice detection preliminary, and spectroscopy or further pho- 
tometry would be worthwhile, but the evidence is certainly as 
strong as for some previous spectroscopic water ice 'detections' 
so we choose to regard this as a confirmed family member for 
the purposes of this paper. This brings the total number of con- 
firmed family members to 10, of the 35 candidate objects. We 
are far more efficient at rejecting candidates though; 8 objects 
have (J - Hs) colours inconsistent with water ice, and cannot 
be true family members. These are 1999 KR16, 1999 RY215, 
2001 QC 298 , 2002 AW 197 , 2003 TH 58 , 2004 PT 107 , 2005 GE 187 
and 2005 UQ513. This is in agreement with Pinilla-Alonso et al. 
[2008], who rejected 2005 UQ513 on the basis of a very red 
slope in an optical spectrum. We also find that 1998 WT31 and 
2001 FU172 have strongly red visible slopes, and can probably 
be rejected as family members without Hawk-I data. Including 
also Makemake and 2004 SB60, which have been shown to lack 
water ice on their surfaces by NIR spectroscopy [Brown et al., 
2007; Schaller & Brown, 2008] and the others listed in the in- 
troduction which have previously been found to have very red 
optical colours, this gives a total of 15 of the 35 candidates that 
are shown not to belong to the family. Finally, we also observed 
1998 HL I5I , 1999 OH 4 , 2000 CG 105 and 2003 HX 56 in the opti- 
cal, but all of these were too faint to put meaningful constraints 
on their family membership. We summarise which objects we 
believe to be family members, which we can rule out, and which 
we do not yet have enough information on in the last column of 
table 4. 
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Table 4. Average colours in BVRIJHs for all candidates (and Eris), and assessment of likely membership based on these colours. 
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Vis. slope 


Ket 


family / 


ft 
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(mag.) 


(mag.) 


(ma: 


§■) 


(mag.) 


(mag.) 


(ve/iUUnrn) 






24835 
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0.65 ± 0.01 


0.39 


±0.01 


0.36 ± 


0.02 


0.65 ± 0.03 


— 


2.0 ± 0.8 


1-8 


Y 




1996 RQ 20 


0.96 ± 0.13 


0.46 


±0.05 


0.71 ± 


0.12 


— 


— 


22.4 ± 6.8 


9,10 


N 


19308 


1996 TO 66 f 


0.68 ± 0.02 


0.39 


±0.01 


0.37 ± 


0.02 


0.61 ± 0.10 


— 


2.9 ± 0.5 


1,10-14 


Y 




1998 HL 151 


0.67 ±0.18 


0.42 


±0.16 


0.79 ± 


0.31 


— 


— 


18.1 ± 16.9 


15,16,32 


9 


181855 


1998 WT 31 


0.76 ± 0.32 


0.51 


±0.25 


0.60 ± 


0.28 


— 


— 


16.6 ± 5.2 


32 


N 




1999 CD 158 


0.83 ± 0.06 


0.51 


±0.05 


0.54 ± 


0.06 


1.38 ± 0.09 


— 


15.8 ± 0.6 


6,32 


N 


40314 


i c\c\c\ is r» i z' 

1999 KR 16 


1 r\i a ao 
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0.75 
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0.74 ± 
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2.99 ± 0.48 
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0.44 ± 
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20.2 ± 35.6 


32 


7 


86047 


1999 OY 3 


A 'H C A Al 

0.75 ± 0.03 


0.26 


± 0.03 


0.33 ± 


0.04 


A OA A 1 O 

0.80 ± 0. 12 


-0.26 ± 0.36 


-0.5 ± 5.3 


3,32 


Y 


86177 


i aaa r> v/ tic 
1999 KY ZD 












A AA i A 1 O 

0.99 ±0.18 


A CO i A TO 

U.jZ + u.zz 




28 


N 




2UUU CCj 105 


1.11 ± 0.25 


0.39 


± 0.13 


0.21 ± 


0.22 






6.7 ± 17.5 


32 


7 




OAA1 T7T T 1 TO 

ZUU1 rU 1/z 


2.91 ± 1.53 


0.66 


± 0.19 


0.53 ± 


0.22 






1A O i ^in A 

39.8 ± 27.9 


32 


N 




onni f\c~ oaq 
ZUU 1 i^C Z^o 


A i A AT 
0.60 ± 0.0/ 
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0.63 ± 


0.07 
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U.jZ + U. 14 
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JN 
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JJJOJ 
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OA T"> n 


IN 




ZUUZ Uh 5Z 


A A 1 i A Hi 

U.91 ± 0.06 


0.66 


± 0.06 


0.56 ± 


0.05 






24.8 ± 4. / 


19,23 


JN 


55636 


2002 TX 300 


f\ s" f- r\ f\r\ 

0.66 ± 0.02 


0.36 


±0.02 


0.32 ± 


0.03 




A A 'I A 1 O 

-0.47 ±0.13 


0.2 ±1.1 


O 1 oo o o 

21,23,32 


Y 


136108 


Haumea' 


l\ /"A A A 1 

0.64 ± 0.01 


0.33 


±0.01 


0.34 ± 


0.01 


A O O . A A 1 

0.88 ± 0.01 


A /'A A 1 1 

-0.60 ± 0.1 1 


-0.6 ± 0.9 


Ol O/l OC ")0 
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Y 




OAA"5 TTV C/^ 
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1 n i 1 it 
1.2/ ± 1.3/ 
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1 O A i 1 H 

lo.4 ± 51. 1 


32 




1ZU1 / 


ZUU j L/r jZ 


u. /u ± u.uj 


0.39 


±0.06 


0.37 ± 


0.05 




AC] _i_ A AO 

-U.D1 ± U.Uo 


1 A -i- 1 1 


v,zo,jz 


v 

1 




700 1" SO 317 












A Ali + 04 


A . A OA 






Y 




2003 TH 58 


0.58 + 0.12 


0.29 


± 0.13 


0.59 ± 


0.15 


-0.13 ± 0.15 


1.29 ± 0.20 


3.6 ± 11.4 


32 


N 




2003 UZ 117 f 














-0.62 ±0.12 


1.1 ± 0.7 rf 


22,27,28,32 


Y 




2004 PT 107 


0.82 + 0.21 


0.65 


± 0.10 


0.68 ± 


0.10 


1.15 ±0.16 


0.54 ± 0.22 


27.9 ± 8.4 


32 


N 




2005 CB 79 


0.73 ± 0.04 


0.37 


± 0.05 


0.36 ± 


0.05 


0.71 ±0.08 


-0.50 ±0.17 


3.1 ±2.5 


32 


Y 


136472 


Makemake 


0.83 ± 0.02 


0.5 : 


+ 0.1" 


0.3 ±o.r 






7.7 ± 8.8 


21,29 


N 




2005 GE 187 












1.22 ±0.19 


0.65 ±0.14 




32 


N 


145453 


2005 RR 43 1 ' 


0.77 ± 0.06 


0.41 


±0.04 


0.29 ± 


0.08 


0.48 ± 0.04 




3.3 ± 6.0 


8,22,26 


Y 


202421 


2005 UQ513 












1.39 ±0.08 


0.30 ±0.12 


18.1 ±2.0 rf 


27,28,30,32 


N 


136199 


Eris 


0.78 ± 0.01 


0.45 


±0.03 


0.33 ± 


0.02 


0.52 ± 0.02 


0.25 ±0.11 


5.9 ±5.6 


21,22,31,32 





Notes: 

a There are no published colours for candidates 1996 TR 66 , 1997 RX,, 1999 OK4, 2000 JG 8 i, 2003 HA 57 . For 2003 QX 91 and 2004 SB 60 we 
measured A-band photometry, but not colours. None of these candidates are included in the table. 

* References: [1] Boehnhardt et al. [2001]; [2] Gil-Hutton & Licandro [2001]; [3] Doressoundiram et al. [2002]; [4] McBride et al. [2003]; [5] 
Tegler & Romanishin [2003]; [6] Delsanti et al. [2004]; [7] Doressoundiram et al. [2007]; [8] Rabinowitz et al. [2008]; [9] lewitt & Luu [2001]; 
[10] Tegler & Romanishin [1998]; [11] Jewitt et al. [1998]; [12] [Barucci et al., 1999]; [13] [Davies et al., 2000]; [14] [Hainaut et al., 2000]; [15] 
[Hainaut & Delsanti, 2002]; [16] [Trujillo & Brown, 2002]; [17] [Sheppard & lewitt, 2002]; [18] Delsanti et al. [2006]; [19] [Santos-Sanz et al., 
2009]; [20] [Doressoundiram et al., 2005]; [21] [Rabinowitz et al., 2007]; [22] [DeMeo et al., 2009]; [23] [Doressoundiram et al., 2005]; [24] 
[Rabinowitz et al., 2006]; [25] [Lacerda et al., 2008]; [26] [Pernaetal., 2009a]; [27] [Pinilla-Alonso et al., 2007] [28] [Alvarez-Candal et al., 
2008] [29] [Ortiz et al., 2007]; [30] [Fornasier et al., 2009] [31] [Brown et al., 2005]; [32] This work. Where colours for a given object are 
published by multiple authors, we quote a weighted mean. 

(R - J) calculated from near simultaneous R and J observations by [17] and [18] respectively. No correction is made to this (or any other 
colour in the table) for possible differences due to changes in rotational phase, although [17] show 1999 KR 16 to have a light-curve amplitude of 
Am = 0. 1 8 mag. 

Although no photometry is published, measurements of the spectral slope for these objects (derived from optical spectra) can be found in the 
table 2 of Fornasier et al. [2009]. We give a weighted mean for each object, and list the references to the original papers in the table. 
" The colours for Makemake are calculated from the BVI photometry from [21] along with the ^-band photometry from [26]. We use a phase 
function of B = 0.05 mag deg~' to correct the iJ-band photometry to zero phase angle, as [21] show that B is approximately constant at this value 
between the V and / bands. The uncertainty is dominated by the uncertainty on the A^band photometry. 
+ These objects listed as confirmed family members by Ragozzine & Brown [2007]. 



Given the high rate of rejection of candidates, we consider 
the likelihood that this is a true family from a statistical point 
of view, or whether the ~ 30% of water ice bodies within 
the candidates could just reflect the proportion within the TN 
region in general. Based on the TNO taxonomy proposed by 
Fulchignoni et al. [2008], the confirmed family members all be- 
long to the BB class, while the rejected candidates come from 
all 4 of the groups (the majority of the newly rejected ones are 
from the red classes RR and IR, since they were mostly rejected 
due to their red slopes). The BB class makes up only 20% of the 



whole TNO population; the proportion of BB within the candi- 
dates (= confirmed family members / candidates « 30%) is high 
but not so unusual given the small numbers of objects involved. 
If instead of taxonomic classes we consider the proportion of 
TNOs with water ice detections (from IR spectroscopy), then 
in the case of the general population we find ~ 50% (from ta- 
ble 1 of Barucci et al. [2008]), making the proportion of bodies 
with water ice in the candidates lower than the general popu- 
lation, although this number contains significant biases as the 
spectroscopy only covers the brightest bodies. We can conclude 
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Fig. 2. Visible and NIR photometry for all candidate family members with observations in at least three bands. The data are nor- 
malized at 0.55^m ( V filter), except in the four cases where no V-band photometry exists. These data are normalised to the /?-band, 
and are indicated by an asterisk next to the designation. The spectrum of Haumea is shown for comparison in each; this is taken 
from Pinilla-Alonso et al. [2009]. The photometry of the large TNO Eris is also shown for comparison; it is not associated with the 
family, and has a spectrum dominated by methane ice [Dumas et al., 2007]. 



that we do not see a significantly larger number of water ice bod- 
ies in the candidate list than in the general population, but this 
ignores grouping in orbital element space. 

In Fig. 4 we show the candidates in terms of their orbital 
parameters semi-major axis, inclination and eccentricity. The 
confirmed family members cluster tightly around the centre of 
the distribution in both plots, where the original orbit of the 
pre-collision Haumea was (Haumea itself now has a higher ec- 
centricity than the centre of the family due to interaction with 
Neptune through orbital resonance [see Ragozzine & Brown, 
2007]). This suggests that the family hypothesis is a valid one, 
but that the spread in orbital elements since the collision is 
less than the range investigated by Ragozzine & Brown [2007]. 
Taking the required collision velocities from that paper ((5v m j n ; 
the minimum ejection velocity required including the effects 
of eccentricity and inclination diffusion in mean-motion reso- 
nances) we find that the largest velocity required by any con- 
firmed family member is 123.3 m/s (for 1995 SM55), while can- 
didates are listed with <5v m ; n up to 250 m/s. If we restrict the 
candidate list to those with Sv m [ n < 150 m/s, we find that the pro- 



portion of confirmed water ice detections rises to 53%, and goes 
up to 64% if we look only at those with lower <Sv m ; n than 1995 
SM55, so the grouping is statistically significant compared with 
the general population of bodies with water ice surfaces within 
the TN region. It should be noted though that there are rejected 
candidates spread evenly across the phase space, including 2005 
UQ513 and 1999 CD158 near to the centre of the family distribu- 
tion, which demonstrates the importance of physical studies of 
the candidates to confirm membership. The remaining unknown 
objects near to the centre of the distribution are 1999 OK4 and 
2003 QX91 (although the latter has high eccentricity and a high 
<5vmin of 222 m/s) which should be high priority targets for fur- 
ther study to measure candidate family member surface proper- 
ties, along with 1997 RX9 which has a low 6v m \ n of 86.8 m/s. 

It is noticeable that the confirmed members remain the larger 
bodies, even though this photometric method is sensitive to wa- 
ter ice absorption on bodies too small for NIR spectroscopy. We 
tested the idea that retention of a water ice surface could be a 
property of only the larger TNOs by looking for a correlation 
between absolute magnitude and the (J - H$) index, but found 
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4.2. Light-curves 



20 40 

vis slope (%/1 00 nm) 

Fig. 3. (J — Hs) colour against visible slope (%/100nm) for 
all candidates (and Eris) and where both measurements have 
been made. Filled circles are confirmed family members, crosses 
show rejected candidates. Haumea itself is the point in the very 
bottom left. 
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Fig. 4. Confirmed family members (grey filled circles), rejected 
candidates (crosses) and those with unknown surface properties 
(open diamonds) plotted in terms of the orbital parameters semi- 
major axis, inclination and eccentricity. Haumea itself is shown 
as a grey circle with a black outline. 



that no such correlation exists. It is likely that there are smaller 
water ice covered family members, however they have yet to be 
discovered or confirmed. We also tested for any correlation of 
the colour with orbital elements and found none; we are dealing 
with a family clustered in dynamical element space, not a conse- 
quence of any correlation of, for example, the presence of water 
ice with semi-major axis. 



This work aimed to test both the membership of the candidate 
family members and also the hypothesis that the family mem- 
bers apart from Haumea itself are composed of almost pure wa- 
ter ice, being made of the reassembled fragments of the outer 
layers of the differentiated proto-Haumea [Brown et al., 2007]. 
This can be probed by testing the density of the family members; 
Haumea is known to have a rock-like density of 2.5-3.3 g cm -3 
(the value found from combining the size from Stansberry et al. 
[2008] and mass from Ragozzine & Brown [2009] agrees with 
the value from the light-curve model of Rabinowitz et al. [2006]) 
but the other family members should have densities at or below 
the density of water ice, ~ 1 g cm -3 . To test this we sought 
to apply the technique of measuring rotation rate and elonga- 
tions using light-curves, which then constrain the density of a 
strengthless body to be 



10.9 a 



(2) 



where a/b is the axial ratio for an ellipsoid and the rotation pe- 
riod P 10t is in hours [Pravec & Harris, 2000]. We can reasonably 
expect the recombined fragments from a collision to be a loosely 
packed 'rubble pile' and therefore strengthless. This method 
only gives a lower limit to the density as the object does not 
need to be spinning at its break up velocity, but must be below 
it, and also that the light-curve amplitude Am only gives a lower 
limit on the elongation, a/b > 10°- 4A " ! . Despite this, when study- 
ing a population a cut off in minimum densities becomes appar- 
ent, which can be used as a reasonable measurement of the bulk 
density of the bodies in the population. This has been clearly 
demonstrated for asteroids [Pravec et al., 2002], where there are 
many (iV > 1000) light-curves available, and also used to derive 
a low bulk density for cometary nuclei [Snodgrass et al., 2006], 
in agreement with the results found by the Deep Impact mission 
[Richardson et al., 2007], despite the relatively low number of 
light-curves available for nuclei. 

We measured partial light-curves for 13 of the candidates 
using EFOSC2, however poor weather during these runs pre- 
vented us from building up the number of light-curves required 
to study the density of these bodies by this statistical technique. 
This was further hampered by the large proportion of the can- 
didates which were eventually rejected as non-family members. 
For 1998 WT31 and 1999 OY3 we have less than 10 points spread 
over three and two nights respectively and there is no significant 
variation. 1999 RY215 and 2003 THgs also show no significant 
variation despite larger data sets. For 1999 CD158, 2000 CG105, 
2002 GH32 and 2005 UQ513 the light-curves show significant 
variation, with ranges of Am = 0.6, 0.45, 0.75 and 0.3 mag., but 
no period could be determined. For 2001 QC298 there is possibly 
a maximum each night in the data, with Am - 0.4, but there can 
be other periods beyond the ~ 12 hours best fit. 2003 SQ317 gives 
a good fit with single peak light-curve of 3.7 hours, while a dou- 
ble peaked light-curve at 7.5 hours also looks reasonable. 2004 
PT107 shows a possible slight variation (0.05 mag), but not a very 
convincing one, with a suggested long period (~ 20 hours). We 
obtained data on 2004 SB60 on two nights which show a varia- 
tion of Am = 0.2, but no clear periodicity. There is a possible so- 
lution at around 17.5 hours, but it is not convincing. 2005 GE187 
has a reasonably convincing single peak light-curve with a 6. 1 
hour period and Am = 0.5. 

The only light-curve in this set of relevance to the density 
of the family members is that of 2003 SQ317, shown in Fig. 5. 
The period of 3.7 hours and the range of Am =1.0 mag im- 
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0.5 1 1.5 2 

rotational phase 

Fig. 5. Light-curve for 2003 SQ317, with data taken on two nights 
(circles from the night of August 30th 2008 and crosses from 
August 31st) folded onto a 3.74 hour period. 

plies a high density, p > 2.0 g cm 4 , however this is for a sin- 
gle peaked light-curve as would be caused by albedo variations 
and not shape. The light-curves of Solar System minor bod- 
ies are more likely to be caused by shape than albedo patterns 
[Jewitt, 2008; Sheppard et al., 2008]. Assuming that this single 
peak light-curve shows half of the period of the true shape con- 
trolled light-curve, the implied density is p > 0.5 g cm 4 , which 
is a weak constraint. Given the sparse light-curve coverage there 
are also other possible periods. We cannot rule out a low density 
and therefore an entirely ice composition for this body. 

Of the other family members 5 of the large bodies with 
confirmed water ice surfaces also have light-curves (not in- 
cluding Haumea itself). 1995 SM55 has a rotation period of 
8.08 hours and Am = 0.19 [Sheppard & Jewitt, 2003] (im- 
plying p > 0.20 g cm 4 ). 1996 TO 66 has P mt = 7.9 hours, 
Am = 0.26 [Sheppard & Jewitt, 2003] (p > 0.22), but both 
the period and amplitude of the light-curve are seen to change 
[Hainautet al., 2000]. 2002 TX300 has a period between 8 and 
12 hours and a low amplitude of Am = 0.08 [Sheppard & Jewitt, 
2003] (p > 0.18). Pernaetal. [2009b] find that 2005 RR43 has 
P rot = 5.08 hours, Am = 0.12 (p > 0.47). Observations from the 
same group find no obvious periodicity for 2003 UZ117. None of 
these light-curves require high densities, although for these very 
large objects it is also likely that the rubble pile assumption will 
be invalidated due to compaction by self gravity, in which case 
finding the density from the light-curve involves assuming fluid 
like behaviour [see Lacerda & Jewitt, 2007]. 

5. Summary 

We have presented optical and/or near infrared colours for 22 of 
the 35 candidate members of Haumea's collisional family that 
were listed by Ragozzine & Brown [2007]. We make use of a 
unique capability of the new Hawk-I instrument at the VLT to 
evaluate the depth of the 1.6 pm water ice absorption band using 
NIR photometry on objects too faint for spectroscopy. We find: 

1. Of the 15 candidates observed with Hawk-I, 7 were found to 
be family members. Most (6) of these were already known 
family members, including Haumea itself, whose confir- 
mation proves the validity of the photometric technique 
used. In addition to the confirmed family members listed by 
Ragozzine & Brown [2007] we confirm the identification by 



Schaller & Brown [2008] of water ice on 2005 CB 79 , and 
identify 2003 SQ317 and as a probable new family member. 
We reject the other 8 candidates observed with Hawk-I as 
interlopers which lack water ice absorption. In general the 
rejected bodies are relatively far from the centre of the family 
in orbital parameter space. 

We present optical colours for 10 candidates and also collect 
all available colour information from the literature for the 
full set. Of the 20 candidates not yet observed with Hawk- 
I there are optical colours for 13. We find that all objects 
where the NIR colour indicates water ice have neutral or blue 
slopes, and consequently we can reject the possibility of wa- 
ter ice on the surface of the very red objects in this sample 
with a reasonable degree of confidence. In this way we rule 
out family membership for a further 5 of the candidates, in 
addition to 2 candidates which are already known to have no 
water ice on their surface from NIR spectroscopy. 

4. Of the 35 family member candidates this gives totals of 10 
confirmed members (29%), 15 non-members (43%) and 10 
that still have to have their surfaces characterised. It appears 
that the family members all fall within the centre of the dy- 
namical region searched by Ragozzine & Brown [2007], so 
we expect that most of the remaining bodies will also be re- 
jected. 

5. We obtained partial /?-band light-curves for 13 of the candi- 
dates, only two of which were subsequently confirmed as a 
family members. Of these 1999 OY3 showed no significant 
variation in the short sequence we were able to obtain on 
it, while 2003 SQ317 shows variations consistent with a 3.74 
hour single peak light-curve, but other periods are possible 
in the sparse data. Neither this nor the existing light-curves 
in the literature for other family members provide strong 
constraints on the density of these bodies, so we cannot yet 
determine whether or not they are 'pure' water ice bodies 
formed from the outer layers of the pre-collision Haumea. 
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